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Abstract 

A known 'sticky-hard-sphere' model, defined starting from a hard-sphere- 
Yukawa potential and taking the limit of infinite amplitude and vanishing 
range with their product remaining constant, is shown to be ill-defined. This 
is because its Hamiltonian (which we call SHS2) leads to an exact second virial 
coefficient which diverges, unlike that of Baxter's original model (SHS1). 

This deficiency has never been observed so far, since the linearization 
implicit in the 'mean spherical approximation' (MSA), within which the model 
is analytically solvable, partly masks such a pathology. To overcome this 
drawback and retain some useful features of SHS2, we propose both a new 
model (SHS3) and a new closure ('modified MSA'), whose combination yields 
an analytic solution formally identical with the SHS2-MSA one. This mapping 
allows to recover many results derived from SHS2, after a re-interpretation 
within a correct framework. Possible developments are finally indicated. 
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I. INTRODUCTION 



In a seminal series of papers [1-3] Baxter first introduced the concept of the so-called 
'sticky hard sphere' (SHS) models, as the simplest - albeit crude - modelization for real 
fluids of spherical particles with a strong surface adhesion. In Baxter's original formulation 
[1-4] and its extension to the multi-component case [5,6] (both hereafter referred to as 
SHS1 model) the pair potential contains - in addition to a hard sphere (HS) repulsion - a 
infinitely deep and narrow attractive square- well, obtained according to a particular limiting 
procedure (Baxter's 'sticky limit') that keeps the second virial coefficient finite [1]. 

Although this model appears rather pathological at first sight, it includes a number 
of interesting features which justify its wide popularity. First, the Ornstein-Zernike (OZ) 
integral equation of the statistical-mechanical theory of fluids can be analytically solved 
for it within the Percus-Yevick (PY) approximation and the solution exhibits a gas-liquid 
transition [1-4]. Second, Baxter's model has already proven to be appropriate for describing 
some properties of colloidal suspensions, micelles, microemulsions and protein solutions with 
short-range interactions as well as some aspects of adsorption, flocculation and percolation 
phenomena, solvent-mediated forces, ionic mixtures, solutions with a small degree of size 
polydispersity and fluids of chain-like molecules (for an illustrative, albeit not exhaustive, 
list of references, see [7-10]). All this means that, in spite of its highly idealized character 
and known shortcomings [7], the SHS1 model is able to capture some important physical 
features of structure, thermodynamics, and phase behaviour of many real systems. 

However, the SHS1 model has its main drawback in its problematic application to mix- 
tures with a large number p of components - as occurs for colloidal suspensions with large 
polydispersity - since this case requires the solution for a set of p(p+ l)/2 coupled quadratic 
equations [5], a task which cannot be accomplished analytically. This important fact has 
originated a more recent attempt of finding an alternative SHS model, which could be an- 
alytically tractable even in the general multi-component case [11-23]. For pure fluids, Brey 
et al. [11] proposed to start from a HS- Yukawa (HSY) potential 
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where r is the distance between particles, (5 = (/cbT) -1 (k B being Boltzmann's constant, and 
T the temperature), a denotes the HS diameter, z the Yukawa inverse range, Eq an energy, 
and T* a reduced temperature (as in SHS1, the factor 12 simplifies subsequent analysis). 
The definition of this second Hamiltonian (SHS2 model in the following) is completed by the 
definition of a 'sticky limit', which in this case amounts to taking z — > +00 [11]. It is worth 
remarking that, unlike its counterpart in the SHS1 model, the starting potential </>HSY( r ) is 
itself independent of temperature. 

For the SHS2 model (in both cases of pure fluids and mixtures), the OZ equation can 
be solved analytically within the mean spherical approximation (MSA) [12,13]. It turns out 
that the SHS2-MSA solution q(r) (or %(r) ) for the Baxter form of the OZ equation has 
exactly the same r-dependence as the SHS1-PY solution, i.e. both are expressed in terms of 
a second-degree polynomial in r. However, the difference between the two solutions lies in 
the density and temperature dependence of their polynomial coefficients (see below), which 
makes the SHS2-MSA solution readily usable even in the multi-component case [14,15,17-19], 
unlike the SHS1-PY one. 

Unfortunately, the common belief that the SHS1 and SHS2 Hamiltonians are different 
but equivalent representations of a unique SHS potential [8,10,14-17,24,25] as well as a 
failure in appreciating the subtle distinction between model and solution [26] have often 
generated a number of misunderstandings and erroneous beliefs in the literature on SHS 
fluids. In particular, both aforesaid PY and MSA solutions have sometimes been regarded 
as corresponding to the same SHS model (Baxter's one) [8,14-17,24]. On the contrary, 
in the present paper we stress that the SHS1-PY and SHS2-MSA solutions are different 
not only as stemming from different closures, but, more importantly, since they refer to 
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different Hamiltonians. In fact, unlike SHS1, the SHS2 model itself is ill-defined from a 
thermodynamical point of view, and thus SHS1 and SHS2 cannot be equivalent. We will 
prove this point by considering the exact second virial coefficient of the HSY potential (1) 
and showing that it diverges in the sticky limit (a very short, preliminary, account of these 
results has been given in Ref. [21]). As we further elaborate below, this pathology is hidden 
in the SHS2-MSA compressibility equation of state (EOS) but is mirrored by a similar 
singular behaviour of the corresponding MSA virial and energy EOS's, never investigated in 
previous studies. 

The second and main goal of the present paper is to propose a new model (SHS3) 
which combines the advantages of both SHS1 and SHS2. As SHS1, it has a finite second 
virial coefficient, and thus it is a well-defined model. As SHS2 however, it admits a simple 
analytical solution within a new closure, referred to as modified mean spherical approximation 
(mMSA) in the following. The remarkable property of the SHS3-mMSA solution for q(r) is 
that it turns out to be formally identical with the SHS2-MSA one, and so are all quantities 
which are immediately derivable from it (such as structural properties and compressibility 
EOS). As a consequence of this mapping, all SHS2-MSA results obtained in the past for 
these quantities can be recovered after an appropriate re- interpret at ion. In addition, we will 
provide new results for other quantities, notably the virial and energy EOS, for which the 
SHS2-MSA solution badly fails. This is not the case in the SHS3-mMSA solution, where the 
energy EOS turns out to be finite, albeit mean-field-like. Finally, although the virial EOS 
displays a singular behaviour within the SHS3-mMSA as well, we will argue that, unlike the 
SHS2 case, this divergence is a consequence of the deficiency of the mMSA closure, and not 
of the SHS3 model itself 

Our findings thus provide a sound theoretical basis for a critical analysis of the existing 
literature on SHS models, and allow to clarify misunderstandings and discard incorrect 
results, previously reported in the literature. 

The plan of the remaining of the paper is the following. In Section II we will briefly 
recall the main features of the SHS2-MSA solution and describe the drawbacks of SHS2 in 
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detail. The new model along with the new closure, which constitute the central part of 
this work, are both discussed in Section III. Final remarks and future perspectives close the 
paper (Section IV). 



II. SHS2 MODEL AND ITS PATHOLOGIES 

For a one-component fluid of N molecules in a volume V, with spherically symmetric 
interactions, the Baxter form of the OZ integral equation [2] is given by 

rc (r) = -q'(r) + 2vrp dtq(t- r) q> (t) , 

(3) 

rh (r) = -q'(r) + 2tt p dt q(t) (r - t) h(\r - t\) , 
where p = N/V denotes the number density, c(r) the direct correlation function (DCF), 
h(r) = g(r) — 1, with g(r) being the radial distribution function (RDF). Moreover, the 
prime denotes differentiation with respect to r. Solving the Baxter equations is tantamount 
to determining the factor correlation function q(r), an auxiliary quantity from which c (r) 
and h (r) can be easily derived. 

An approximate integral equation can be obtained by adding to the OZ equation some 
approximate 'closure' relating c(r), h (r) and the potential 0(r). 



A. MSA solution 

For potentials with a hard-core part, the mean spherical approximation (MSA) reads 

|-[l +7 ( r )], 0<r<a, 
cmsa (r) = <^ (4) 
( -/30taii (r) , r > a, 

where 7 (r) = p J dr' c (r') h (|r — r'|), tail (r) is the potential outside the core. Since g (r) = 

1 + 7 (r) + c (r) , the MSA may also be written as 

fl-MSA (r) = e HS (r) [1 + 7 (r) - /30 ta ii (?")] » ( 5 ) 

where ens ( r ) = ex P [~ /30Hs( r )] = & ij' — °~) , with 0hs(^) being the HS potential, and 9(x) 
the Heaviside step function (9(x) = for x < 0, and 6*(a;) = 1 for x > 0). 
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For the HSY potential the MSA closure becomes 

CHSY-MSA(r) = K e-< r ~^ /r, r > a, (6) 

where the parameters are the same as in Eq. (1). The analytic MSA solution gHSY-MSAfj), 
for the HSY fluid [27] yields in the sticky limit (z -> +oo, ) the SHS2-MSA solution 

1 



gSHS2-MSA(V) 



^a(r 2 - a 2 ) + ba(r - a) + q a & 2 



8(a-r), for r > 0, (7) 
1 + 2i] 12q a rj , 3?7 6^77 



(8) 
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q a = K /a 2 = e* . (9) 

where r] = (n/6)pa 3 is the packing fraction. As remarked, the expression (7) is formally 
identical with the SHS1-PY solution with the crucial difference that gSHS2-MSA dep^,-^ on \y 
on temperature (being proportional to T* _1 ), whereas g^Hsi-PY j^gj^g on temperature 
and density (in a more complex way) [2]. 

Once q{r) is known, all structural and thermodynamic properties can, in principle, be 
calculated. It is worth mentioning that, with a few exceptions [12,16,20,22,23,25], the ther- 
modynamic properties of SHS2 are still mostly unexplored. However, we stress once again 
that the main aim of our paper is not to present new results on the SHS2-MSA thermody- 
namics (although this task will be accomplished too), but to reveal a dramatic fault of the 
SHS2 potential itself (irrespectively of any approximate closure), which, surprisingly enough, 
has never been observed in the previous literature, but will emerge from the following simple 
analysis. 



B. Sticky limit of the exact second virial coefficient of the HSY fluid 

The Hamiltonian of the SHS1 model was introduced by Baxter as a limiting case, through 
a clever definition of both a starting potential and a limit procedure in such a way that the 
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contribution of the vanishing square- well tail to the second virial coeffient, B 2 , remains 
finite and nonzero. It is instructive to consider the behaviour of B 2 for the SHS2 model, by 
considering the result of the z — > +00 limit for the exact second virial coefficient _gHSY-cxact 
of the HSY fluid. This can be calculated from the general definition, 

B 2 = -2tt / dr r 2 f{r), (10) 
J 

with /(r) = e(r) — 1, and e(r) = exp [—(3(f) (r)] being the Boltzmann factor. One finds 

^HSY-cxact _ ^[^3 _|_ ^^HSY-cxact 

3 

where 

POO * r -. ^ 

Aj gHSY-cxact = _ 2% J dr r 2 j exp e ~z{r-c) ^ _ j j ^ 

It is now easy to show that A^HSY-cxact diverges j n ^he z — > +00 limit. To this aim we note 
that, for x = zK e'^'^/r > 0, one has e x — 1 > x + x 2 /2 and hence we can use the bound 

jH dr r 2 ( e x - 1) > jH rfr r 2 + ^ = K (<r + ^) + ^ (11) 

As the right hand side of Eq.(ll) diverges as 2 — > +00, we have thus shown that 

R S HS2-exact = ^ pHSY-exact = _ QO> 

This result is exact and independent of any closure, and reflects an inconsistency of the HSY 
potential with the definition of sticky limit employed when setting up the SHS2 Hamiltonian 
[11,12]. As a consequence, the SHS2 Hamiltonian is ill-defined from the outset, and the 
corresponding model (which cannot be a different representation of Baxter's one) must be 
discarded. 

This is, however, a surprising result in some respect. One may rightfully wonder why 
no trace of the pathological nature of the SHS2 model has ever been revealed by a number 
of structural studies [14,15,17-19] carried out on its multi-component version. As discussed 
in Section IV, all these structural results based upon the one and multi-component SHS2- 
MSA solutions are, in fact, fully correct after their re-interpretation in terms of SHS3. Before 



doing this, however, it is instructive to consider those SHS2-MSA thermodynamic properties 
of the one-component fluid which have not been investigated so far. This will display the 
pathological nature of the model as shown next. 



C. Sticky limit of the MSA equations of state for the HSY fluid 



Mier-y-Teran et al. [12] obtained, for Z = (3P/p, the MSA compressibility (C) EOS 
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by integrating with respect to density the compressibility equation [2,28] 

(fV [1 ~ 2ip?(0)|W ' (14) 

where q(k) denotes the unidimensional Fourier transform of q(r). 

On the other hand, no expressions for the SHS2-MSA virial and energy EOS were given 
either in Ref. [12] or in subsequent literature on this subject, to the best of our knowledge. 
We thus tackle this analysis here. 

The most direct way we have followed is to consider the z — > +00 limit of the HSY- 
MSA virial and energy pressures. Convenient expressions for these quantities were given 
by Cummings and Smith (Eqs. (18)-(20) of Ref. [29]). Unfortunately, their sticky limit 
requires a rather elaborate analysis, which we will not report here. However, the final result 
is rather simple: both the virial and the energy MSA pressures of the HSY fluid diverge in 
the sticky limit. For the MSA energy EOS, the same conclusion can be drawn by exploiting 
the alternative form reported by Herrera et al. [30]. It is worth stressing again that this 
important feature has never been pointed out before, to our knowledge. 

To convince the reader that our statement is correct, we now provide a simpler demon- 
stration, based upon the analysis of the sticky limit of the MSA second virial coefficients of 
the HSY fluid, as obtained from the density expansion of ^ HSY-MSA ) ^ ^hsy-msa^ anci 
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^hsy-msa^ j g we jj known that, in any approximate theory, compressibility, virial and 
energy EOS's may yield different virial coefficients, and the first one at which this difference 
begins to appear depends on the chosen closure. We now show that, in the SHS2-MSA case, 
a non-consistency already appears at the f?2- level, and partly veils the singular character 



of the SHS2 Hamiltonian. Nonetheless, as it will be shown below, both (b 



hsy-msa\ 



IV 



and 



(^Bf SY ~ MSA ^j are divergent, and this strongly supports the aforementioned statement about 
the singular character of the virial and energy EOS for the SHS2-MSA solution. 

The quantity ^B^ SY_MSA ) can be easily computed directly from the density expansion 
of d(3P/dp, given by Eq.(14), upon using Eqs. (8)- (9). To evaluate (s^ - ^) and 
(-S^ SY_MSA ) , it proves convenient to exploit the following expressions for the MSA virial 
and energy EOS of a HSY fluid [31,32] 



L V 



1 + 4?7 3hsy(^ + ) -I 



(15) 



Z| SY = Z HS + A V \ - g^ Y (a + ) - fl&(<7- 



-I 



roo 

I = K a' 3 dr rg mY (r) (1 + zr) e'^'^. 

J a 



(16) 



(17) 



Using the low-density expansion of the MSA approximation to <7hsy(?")> Eq. (5), one gets 

1 \ 1 
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(19) 
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where Vq = (ir/6)a 3 is the particle volume. The difference among the three results for 
^hsy-msa j g a c j ear manifestation of the known thermodynamic inconsistency of the MSA 
[29,33]. It is however quite surprising that discrepancies already appear at the £> 2 -level. 
Note that (S 2 HSY - MSA ) C > (B* SY ~ MSA ) E > y „ 
with increasing z, and become dramatic in the sticky limit which yields 



j^hsy-msa^ _ Thege differences magni f y 
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(£f S2 - MSA ) c = ^(l-^), (21) 

^SHS2-MSA^ =.00= ^SHS2-MSA^ (22) 

Eq. (22) thus confirms our previous statement about the divergence of the MSA virial and 
energy pressures of the HSY fluid in the sticky limit. On the other hand, Eq. (21) gives 
a finite value for ^|hs2-msa^ j n a g reemen ^ w ith the corresponding compressibility EOS. 
However, we stress that this result is due to the fact that the MSA involves an approximation 
at the level of Boltzmann factor [34], i.e. 

e MS A (r) = 6(r - <r) [1 - /30 tail (r)] , (23) 

equivalent to the linearization exp [— /90 ta ii (?")] — 1 — P<Pta,n (r) (as can be inferred from a 
comparison between the zeroth-order terms in the density expansion of g cxa ct( r ) and (?msa ( r ), 
i-e- 5-exLt ( r ) = e(r) = e HS (^) exp [-/30 tail (r)] and g$ SA (r) = e HS (^) [1 - /3<Ataii respec- 
tively). The replacement of the exact e H sv(^) in the expression for -B^ SY with its MSA 
counterpart (23) just leads to the ^g| HS2 ~ MSA ^ result, in the sticky limit. 

Summarizing, we can conclude that the existence of the MSA solution, <7shs2-msa(?")) is a 
fortuitous consequence of the fact that the divergent character of the SHS2 model is masked 
by the MSA linearization (23). 

D. MSA internal energy 

As a final point of this first part, we wish to point out an interesting different behaviour 
between the internal energy in the SHS2-MSA and SHS1-PY solution. 
Let us now consider the energy (E) equation, 

u cx = U cx /N = 2vrp / air r 2 <?(r)0(r), (24) 
Jo 

where U ex is the excess internal energy. Within the MSA, one gets for our HSY fluid 
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-/3< SY - MSA = 27T p / dr r 2 [1 + thsy-msa (r) - (3cf> Y (r)] [-/tyy (r)] 

{/*oo 
| rfr r [1 + 7hsy-msa (r)] -ze^^ + -^z 

Since 7hsy-msa(?") is finite and sufficiently regular everywhere, it is easy to show, integrating 
by parts or using the property lim 2 ^ +00 ze~ z ( r ~ a ^9(r — a) = 5 + (r — a) [35], that the last 
integral remains finite as z — > +oo. In conclusion, -u^f Y_MSA diverges linearly. This should 
be compared with the logarithmic divergence found in SHS1-PY energy [1]. We speculate 
that this difference in the internal energy behaviour parallels that in the corresponding 
EOS. It is then clear that, in the SHS1-PY case, such a weak divergence does not in itself 
constitute an impediment to the existence of finite EOS's in the sticky limit [1,7,36]. 

As a side comment, it is worth noticing that our finding contradicts Ref. [37], where an 
MSA energy with a non-singular sticky contribution (Eq. (42) of that paper) is reported . 

III. NEW MODEL AND NEW CLOSURE 

Given the above premise, one could, at this point, suspect that all previous findings 
based upon the SHS2 model (and notably, those referring to polydisperse colloidal fluids 
[18-20]) should be discarded. This is not so, and all those results are, in fact, correct. The 
second aim of our paper is to show that an analytic solution with the same functional form as 
<7shs2-msa(?") is obtained within a new and well-defined model (which we call SHS3), coupled 
with a simple new closure. As a consequence, most of the results derived from SHS2 can be 
recovered after an appropriate re-interpretation. 

A. SHS3 potential 

We define a new SHS Hamiltonian, which represents the simplest correct alternative to 
SHS2, and is analytically solvable within a novel closure. The basic idea hinges on Baxter's 
trick of a logarithmic tail, combined with the infinite-ranged Yukawa expression of SHS2. 
The starting potential (hereafter referred to as M3) is 
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+oo, < r < a, 

P<pMs(r) = { (25) 

- In [l + zK e-'^''-^ /r], r> a, 

where K depends on T according to Eq. (2), and the sticky limit corresponds again to 
z — > +oo. 

We note that at large r values, the quantity x = — f3(f) Y (r) = zK e _2<r ~°")/ r > appearing 
in the argument of the logarithm, becomes so small that the approximation ln(l + x) ~ x 
can be used. In other words, M3 has the same large-r asymptotic behaviour as — /30y (r), 
but differs from the Yukawa tail near contact. This is depicted in Figure 1, where it can be 
clearly seen that this difference between the two potentials increases as z increases. However, 
for the M3 model, the Boltzmann factor reads 

e M3 (r) = 9 (r - a) [l + zK e~ z ^/r] , (26) 

with the essential consequence that the corresponding exact second virial coefficient remains 
finite in the sticky limit, i.e., 



^MS-exact = ^ 



1 - 5 ■ * ( 1 + — 

za 



(27) 



^SHSS-cxact = j^^MS-exact = ^ ^ _ _ ^ 

Note that this ^-expression, exact for SHS3, is formally identical with £§hs2-msa which 
is only approximate for SHS2. 

As a further check that we are on the right track, we have also computed the exact SHS3 
third virial coefficient 

RSHS3-exact = ^MS-exact = 2 ( w _ A + J_ 1 ), ( 2 9) 

3 z^+oo 3 V T* T* 2 18T* 3 / 

It is remarkable that #| HS3 - cxact coincides with £| HS1 ~ cxact ; once t* is replaced with its 
counterpart r of Baxter's model [1]. 
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B. Modified MSA closure and solution 



We define a modified mean spherical approximation (mMSA) 



CmMSA (r) = < 



/tan (r) 



-[1 + 7(0], 

exp [-/30 tail (r)] - 1, 



< r < a, 



r > a. 



(30) 



In terms of RDF, this closure reads 



9mMSA (r) 



e HS (r) [1+7 (r) + /tail (r)] 



e (r) + e HS (r) 7 (0 , 



(31) 



since e H s (0 [1 + /tail (r)] = e (r). 

Although this closure is not completely new (see Ref.s [38,39]), it was never formulated 
in the present general form. Yet it is a very natural choice, since it yields the correct 
zeroth-order term in the density expansion of the DCF: c exa ct(^) — ► fif) when p — > 0. 

The advantage of coupling the M3 potential with the mMSA is that one finds 
CM3-mMSA(?") = / (r) = zKq e~ z ( r ~ a )jr (r > a), which is identical with the usual Yukawa 
closure. This mapping allows the immediate identification between the solution gM3-mMSA( r ) 
with (/hsy-msa^), an d thus in the limit z — > +oo we get 



In short, the solution given by Eqs. (7)-(9) is recovered, but within a new well-defined model, 
with a finite second virial coefficient. 



<?SHS3-mMSA(?") = 9SHS2-MSa(^)- 



(32) 



C. Equations of state 



1. Compressibility route 



Since Z^ HS3 mMSA can be computed directly from <ZsHS3-mMSA(?") an d Eq. (14) (in fact, 
sticky limit and p-integration commutate), it results that 



13 



^•SHS3-mMSA _ ^SHS2— MSA 



SHS3-exact 
2 > 



Expanding this quantity in powers of p, one finds that (_Bf HS3 mMSA ^ — B 
whereas (flfHS3-mMSA) = y 2 [ 1Q _ ^ ( T *yi + i (j>*y^ differs from £|HS3-e*act . 

2. Virial route 

Let us insert Eq. (31) into the virial (V) equation, 

O77- roo 2"7r 
Z v = l + —p dr r 3 g(r) [-/ty'(r)l = 1 + —pJy (34) 
3 jo 3 

The integral can be written as Jy = J± + J 2 , with 

Jx = / dr r 3 e(r) [— /30'(r)l = / dr r 3 e'(r) = / dr r 3 /"'(r), 

JO JO JO 

■h = \ dr e HS (r) r 3 7 (r) [-/30'(r)] 

JO 

e HS (r) r 3 7 (r) [-/¥hs(0 - /tyUr)] . 





Integrating Jx by parts, and observing that the boundary terms vanish, one finds that 
(27r/3) Ji = B 2 . On the other hand, if J 2 is written as J 2a + J 2 fe, with J 2a = 

j* °°dr e H s (r) r 3 7 (r) [-/30ns( r )L and usin S e ns (r) h/30Hs( r )] = e hs ( r ) = <K r - °0> one 
gets 

27r r z* 00 1 

z mMSA = 1 + B2f) + _ p | a 3 7mMsA + I rfr (r) [-/ty^)] | . ( 35 ) 

In the particular case of the M3 potential, — f34>' ta x\{r) assumes increasingly large negative 
values as z increases, whereas 7M3-mMSA( r ; z ) < 7M3-mMSA(c; z ) remains bounded, in such 
a way that the integral J 2 j, diverges in the sticky limit (for shortness, our detailed analysis 
of this point is not reported here). 

Eq. (35) and these results imply that the mMSA-virial EOS diverges. However, since the 
second virial coefficient is finite and exact, and (#SHS3-mMSA^ diverges whereas ^Hss-cxact 
is finite, the singularity of ^s. HS3 - mMSA 1S no t due to the Hamiltonian, but is surely caused 
by the mMSA closure. 
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3. Energy route 



Finally let us explore the energy (or free energy) route for temperature-dependent po- 
tentials. It can be shown [3] that, for any potential, 



d_ 



P(A-A id ) 
N 



9 2 , v d \P4> (r)] 

9 Jo r g ^ — d( — 



(36) 



where ( denotes an arbitrary parameter from which (3<f> (r) depends upon, and A- ld indicates 
the ideal gas free energy. In particular, when ( = (3 and g{r) = g m MSA ( r ), one gets 



d_ 

dp 



P(A-A id ) 
N 



= 2np 



mMSA \ " u ^ 

For the M3 potential, this expression becomes 

P (^M3-mMSA — A id ) 



deir) r° , ,, 
dr r — r^- + / dr r 7 mM SA (r) 



(9 



tail 



(r)]" 
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N 



-12 V e 



j J_ 7M3-mMSA /3, g) 



+ (V 2 ) 



, (37) 



With the boundary condition T — > oo (/3 — > 0), corresponding to the HS case, this equation 
can be integrated with respect to and in the z — > +oo limit one finds 



/3(A SH S3-mMSA - A HS )/N = ~J^V, 



(38) 



(since sticky limit and /^-integration commutate). From Z = rjd(pA/N)/dr] we then get 



7 SHS3-mMSA 



= Z HS - 



(39) 



It is noteworthy that these results for the Helmholtz free energy and the energy-EOS are 
van der Waals-like. The effect of the surface adhesion enters these expressions only at the 
level of the second virial coefficient (which is exact), while all higher-order virial coefficients 
predicted by Eq. (39) coincide with the pure HS ones. The reason for this may be traced 
back to the fact that the mMSA closure takes into account only the zeroth-order term in 
the density expansion of the DCF outside the core. 
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IV. CONCLUSIVE REMARKS 



Three different models of 'sticky hard spheres' have been treated in this paper. The first 
two (SHS1 and SHS2) were already present in the literature, while the last one (SHS3) is 
new and constitutes the main contribution of the present work. 

All these models describe a fluid of rigid spherical particles with infinitely strong surface 
adhesion, defined through an appropriate 'sticky limit' which constitutes an essential part of 
the model. The choice of the starting potential strongly influences that of the approximate 
closure to be used to solve the OZ equation analytically. In SHS1 a suitably defined square- 
well allows the solution within the PY approximation, whereas in SHS2 the starting point is 
a HS- Yukawa potential which requires the MSA approximation. However, the most crucial 
point is that, while SHS1 is a perfectly well-defined model, SHS2 is not, since its exact second 
virial coefficient is divergent, as we have proven in this paper. This singularity is mirrored 
by a similar behaviour appearing in both the virial and energy SHS2-MSA EOS's, but it 
should be emphasized that it is a weakness of the model itself and is not due to the MSA 
approximation. This means that the SHS2 model itself is not properly defined and must be 
abandoned. 

To replace it, a proposal has been put forward in the second part of the paper. Our 
recipe is based upon the introduction of a new potential (M3), coupled with a new closure 
(mMSA), in such a way that the Baxter-OZ equations turn out to be analytically solvable, 
with formally the same solution q(r) found in the SHS2-MSA case. This remarkable corre- 
spondence allows a recovery of all previous results for structure factors and compressibility 
EOS based upon the SHS2-MSA solution [14,18-20], after a re-interpretation in terms of 
SHS3-mMSA. 

The SHS3 model is the simplest correct alternative to SHS2. It is well-defined, since the 
sticky limit of the exact second virial coefficient, £?SHS3-exact^ rema j ns fi n n e . Furthermore, 
its starting potential, M3, is asymptotically equivalent to the HSY one at large r values, 
being different from it only in the contact region. 
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As regards the closure, the mMSA is correct in the zero-density limit, and has, conse- 
quently, a higher thermodynamic consistency than the MSA. As a matter of fact, compress- 
ibility, virial and energy routes all generate the same, exact, B 2 . On the other hand, one 
cannot expect from the mMSA to go further in consistency, since this closure is still rather 
poor, taking into account only the zeroth-order term in the density expansion of the DCF 
outside the core. We have shown that discrepancies are already found at the level of third 
virial coefficient: (^fHss-mMSA^ is temperature-dependent but different from sfHSs-exact^ 
j^SHSs-mMSA^ diverges, while ^SHSs-mMSA^ ig temperature-independent and equal to its 
HS counterpart. Clearly, these three results could coincide within a more refined (density- 
dependent) closure. 

It would be very interesting to investigate whether the SHS3 model admits further an- 
alytic solutions, within more sophisticated approximations (such as, in particular, the PY 
approximation). We plan to do such analysis in future work, including also a more detailed 
comparison with Baxter's original model. 

Finally, it is worth noting that the extension of the SHS3 model to mixtures can be easily 
carried out. We hope that the corresponding solution can provide a simple useful tool for 
further studies on structural and thermodynamic properties of polydisperse colloidal fluids. 
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FIG. 1. Comparison between /34>M3(r) (solid curves) and /30hsy(j) (dashed curves) at 
T* = 0.17, for increasing values of the dimensionless inverse-range parameter, z* = za. Note 
the difference in both vertical and horizontal scales between the two parts of the figure. 
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